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Abstract 

Q-f We develop a model for 7-ray emission from the outer magnetosphere of pul- 

O ' sars (the outer-gap model). The charge depletion causes a large electric field 

■ which accelerates electrons and positrons. We solve the electric field with ra- 

diation and pair creation processes self-consistently, and calculate curvature 
>- ! spectrum and Inverse-Compton (IC) spectrum. We apply this theory to PSR 

B0833-45 (Vela) and B1706-44 for which their surface magnetic fields, ob- 
served thermal X-rays are similar to each other. We find that each observed 
cut-off energies of the 7-rays are well explained. By inclusion of emission 
outside the gap, the spectrum is in better agreement with the observations 
than the spectrum arising only from the inside of the gap. The expected 
TeV fluxes are much smaller than that observed by CANGAROO group in 
the direction of B 1706-44. 

1 Introduction 

The Compton Gamma Ray Observatory has detected seven 7-ray pulsars. 
The observed light curves and energy spectra have been used to discriminate 
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possible radiation models. Moreover, the next-generation 7-ray space tele- 
scopes and the ground based Cherenkov telescopes will further constrain to 
the models. 

The pulsed 7-rays imply that the electrons and positrons are accelerated 
up to about 10 TeV in the magnetosphere. For the outer-gap model (Chen, 
Ho & Ruderman, 1986), Hirotani & Shibata (1999, HS) solved the accel- 
erating electric field with curvature radiation and pair creation processes 
self-consistently. They showed that the electric field along the magnetic field 
does not extend to the light cylinder, and calculated curvature spectra can 
explain the EGRET observations in the GeV bands. However, the gap emis- 
sion could not explain the observations in the MeV band. Hence, to improve 
the HS model, we take account of the radiation from the outside of the gap 
and compare the corrected spectrum with the observations. 

2 One dimensional model 

In this section, we introduce the HS model and represent the Poisson eq. 
which describes the accelerating electric field, the continuity eqs. for particles 
and 7-rays , following HS. 

2.1 Basic equations 

We deal with the structure along the magnetic field lines. If the gap width 
(W) along the magnetic field is much less than the light radius (wi c ), we 
can approximate the magnetic field lines as straight lines in the gap, and 
the electric field structures can be treated as one-dimensional, where the arc 
length from the surface along the last-closed line is denoted by s : we can 
write down the Poisson equation as 

^i= 47re fjv + -jv_-^V (i) 

us V e J 

where En is the electric field along the magnetic field, N+ (AL) is the elec- 
tron (positron) number density, and pcj is the Goldreich- Julian (GJ) charge 
density. Above equation describes that the charge depletion relative to pcj 
causes E\\. 
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By taking account of the electron-positron pair creation process, the con- 
tinuity equations for particles and 7-rays are 



±b4-< n * 



roc 

/ da r [ri P G + + ri-G-], (2) 
Jo 



ds V B ) c cos \l/ Jo 

d (G±\ -y P ±G± + r] c N± 

±b Ts\-b) = ^ol* ' (3) 

where G + (GJ) is the distribution function of outward (inward) propagating 
7-rays, e refers to the photon energy in units of the electron's rest mass 
energy, B is the magnetic field strength, i] p ± is the pair-creation rate, r\ c is 
the emissivity of curvature radiation, and ^ is the angle between the particle's 
motion and the meridional plane. We describe G± in several energy bins and 
represent then as G± (i = 1,2, ...). The particle's Lorentz factor in the gap 
is obtained by assuming that particle's motion immediately saturated at the 
balance between the electric and the radiation reaction forces, i.e., 



sat 



2e 



, (4) 



where R c is the curvature radius of dipole magnetic field lines. 

We impose some boundary conditions. The inner (si) and outer (^2) 
boundaries are defined so that E\\ vanishes, i.e., E\\(si) = £j|(s2) = 0. We 
assume that the 7-rays do not come into the gap through the boundaries, 
G + (si) = G l _(s2) — (i — 1,2, ....m). We allow the particles to come into 
the gap. The particle flux is given by the non-dimensional parameters ji and 
j 2 as follows: 

N + ( Sl ) N.(s 2 ) 

■71' r>p/„ \ J7T— =J2- (5) 



nB( Sl )/27rce J ' ns(s 2 )/27rce 
The particle continuity equation (Q) yields 

Njs) N_(s) 

+ , w = const along s = j tot . (o) 



nB(s)/2nce VtB(s)/2Tice 

From eqs. (|3p and (^), the current carriers created in the gap per unit flux 
tube is jg a p = jtot — ji — 32- The total current should be determined by the 
global condition which includes pulsar wind. Therefore, we use (jtot, ji, 32) 
as free model parameters. 
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Table 1: Observed Parameters 



pulsar 


distance 




B\2 


kT s 


Ref. 




kpc 


rad/s 


1Q 12 G 


eV 




Vela 


0.5 


70.6 


3.4 


150 


Ogelman et. al. 


B 1706-44 


1.8(DM)/2.5(HI) 


61.6 


3.1 


143 


Gotthelf et.al. 



DM : Dispersion Measure , HI : HI absorption 



2.2 X-ray & infrared (IR) field 

Because the 77 pair-creation process is important in the outer magneto- 
sphere, we need the X-rays for the target-photons in our case. In the present 
paper, we use the observed black body radiation from the pulsar surface for 
the X-rays (Table 1). We also need the IR field to calculate IC flux. The 
IR field is inferred from optical and radio observations for Vela, X-ray and 
radio observations for B1706-44 with a single power-low, because there are 
no available IR observations. 

3 Radiation from outside of the gap 

Near the boundaries, the real Lorentz factor of the accelerated particles must 
be lager than T sat given by eq.(|J), because the particle's cooling time for the 
radiation becomes larger than the crossing time for the gap. So, the particles 
come out from the gap with T out ~ a few times 10 7 and emit 7-rays outside 
of the gap. Since the typical damping length for the curvature radiation is 

Sm e c 2 R 2 c _ f n vVry 3 / r c \ 2 

loam - 2^iT - °" 4 ^ Uorads-i J [WJ UtoJ ' (?) 

the 7-ray radiation from the outside of the gap are also important unless 
W ~ w ic . 

We apply the our model to Vela and Bl 706-44 with the observed param- 
eters in Table 1. We adopt (j tot , j u j 2 ) = (0.201,0.191,0.001) and a inc = 45° 
for the angle between axes of rotation and magnetization for both pulsars. 



4 



8000 




distance along field line (10 A 6 m) 

Figure 1: The accelerating field for Vela (solid-line) and Bl 706-44 (dotted- 
line for dis=1.8kpc and dashed-line for dis=2.5kpc). (jtot, ji, j-i) — 
(0.201,0.191,0.001) and a inc =45 deg. 

4 Result 

4.1 Electric field structure 

The calculated En for both pulsars are shown in Fig.l. The gap width W 
is shorter than w\ c . W is characterized by the pair-creation mean free path, 
which is given approximately by Wocc/ (J r q p ^r] c de) 1 ^ 2 by using the fact r) p+ <^ 
rj p -, due to the difference in the collision angle between 7-ray and X-ray. For 
Vela and B1706-44, one finds the mean free path to be shorter than w\ c . 

We find that Vela has a nearer distance from the surface to the gap and 
larger calculated E\\ than B1706-44 when we adopt the same {jtot, ji, jz) an d 
cii nc . This is because Vela has the shorter rotation period and larger GJ 
density in the gap than B 1706-44. The dependence of the assumed distance 
from the earth to B 1706-44 is also shown in Fig.l. In general, if we adopt 
a nearer distance to the pulsar, the electric field becomes large, because the 
decrease in the estimated X-ray luminosity from the observations extends the 
gap width. 
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Figure 2: (a):7-ray spectrum for Vela. The total spectrum (solid-line) in- 
cludes the radiation from the outside of the gap (dots-line) as well as gap 
emission (dashed-line). (b):Total 7-ray spectrum for B1706-44. The depen- 
dence of distance is shown. The IC spectrum is also shown in figure. 



4.2 Gamma-ray spectra 

The calculated spectrum of outward propagating 7-rays radiated in the gap 
for Vela is shown in Fig. 2 (a) as dashed-line. The spectral cut-off around 
GeV is responsible for the acceleration limit. We find that this spectrum 
is in agreement with the EGRET observations (Thompson et.al. 1999) in 
the GeV bands. In the MeV bands, however, it is inconsistent with the 
observations. This is because the value of T sat makes the curvature spectrum 
with E 2 F oc E a ,a ~ 4/3 in the MeV bands, although the observations have 
a ~ 1/3. 

In §3, we have pointed out that if W -C ca; C) the curvature radiation from 
the outside of the gap is important. In the outside of the gap where E\\ is 
vanished, the particles lose their energy by the radiation, the spectrum of 
which extends to the MeV band (dotted-line in Fig. 2 (a)). By inclusion of 
this emission, the total spectrum (solid-line in Fig. 2 (a)) is in good agreement 
with the EGRET observations. 

The calculated total 7-ray spectrum for Bl 706-44 is shown in Fig. 2 (b). 
We find that the calculated peak energy becomes slightly less than Vela 
because E^ 1706 < E^ ela (§§4.1), and this peak energy also explains the ob- 
servations. As mentioned in §§4.1, since the calculated E\\ becomes large as 
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we assume the nearer distance to the pulsar, the spectrum becomes hard if 
we adopt a nearer distance. For 1.8 kpc, the calculated spectrum appears 
to be consistent with the observations. However, we must assume a large 
cross-section area of the gap, A± = (10W) 2 ~ (ro/ c ) 2 , to obtain the observed 
fluxes. Because the gap locates at s ~ 0.5zui c (Fig.l), such A± should be 
unrealistic. 

Fig. 2 (b) also shows the calculated IC spectrum from the gap. The sharp- 
cut off in the spectrum which corresponds to the acceleration limit appears 
around lOTeV. Since IR flux might be overestimated, TeV flux would be less 
than ~ 10~ 12 erg cm~ 2 s" 1 . Moreover, this calculated flux hardly depends on 
the distance to the pulsar. On these ground, we conclude that it is difficult to 
explain the unidentified TeV components observed by CANGAROO group 
(Kushida et.al., 2002) in the direction of B1706-44 with this model. 

5 Discussion 

In summary, we obtained the spectrum in good agreement with the EGRET 
observations for Vela by inclusion of the curvature radiation from the outside 
of the gap. We found that the observed peak energies for Vela and B1706-44 
may imply that the almost the same currents in units of GJ value are running 
through the gap for both pulsars. 

From Fig.2 (a), we recognize that the calculated spectrum is inconsistent 
with the COMPTEL observations. If we try to explain this observations, we 
need the very small curvature radius as compared with the dipole, which is 
unlikely. Therefore, this MeV emission will be obtained by inclusion of the 
synchrotron emission by pairs. 

In §§4.2, we showed that we need very large A±_ to explain the observations 
of B1706-44. This may be due to the small jtot, about 20% of the GJ current. 
However, if we adopt the nearly (j tot ~ 1) or super (J tat > 1) GJ current with 
this model, the gap width will be quenched. 

Quite recently, Hirotani et.al. (2002) showed that the value of the T sat 
given by eq.(^) and also 7-ray fluxes calculated with this T sat are overesti- 
mated, because particles in the gap do not immediately saturate. But, the 
general features of the gap model and the result that we must take account 
of the radiation from the outside of the gap to explain the observations above 
lOOMeV are not altered. 
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